Nondemolition measurements of a single quantum spin using Josephson oscillations. 
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We consider a Josephson junction containing a single localized spin 1/2 between conventional 
singlet superconducting electrodes. We study the spin dynamics and measurements when a dc- 
magnetic field B || z acts on the spin and the junction is embedded into a dissipative circuit. We 
show that when tunneling or a voltage are turned on at time t = the Josephson current starts 
to oscillate with an amplitude depending on the initial {t — 0) value of the spin z-component, 
Sz = ±1/2. At low temperatures, when effects of quasiparticles may be neglected, this procedure 
realizes a quantum-non-demolition (QND) measurement of Sz. 

PACS numbers: 74.50. -Hr, 03.67.-a, 03.65.Yz 

Quantum measurements in mesoscopic systems by use 
of tunneling attracted recently great interest. One of 
the motivations is the challenge of the single spin detec- 
tion [1-3]. The other motivation comes from the quan- 
tum computing where the final state (after the computa- 
tion) of a qubit must be measured. The meter couples, 
usually, to a single observable of the qubit, e.g., az- If 
this observable commutes with the spin Hamiltonian, the 
QND regime is realized, i.e., the two possible eigenvalues 
of az can be measured with proper probabilities even by 
a weakly coupled meter. Otherwise the initial state is 
quickly destroyed and one can only observe the steady 
state properties of the qubit and the meter performing 
continuous measurements [4-6]. In the spin detecting 
tunneling schemes all components of the localized spin 
are coupled to the tunneling electrons via the exchange 
interaction. Thus QND measurements seem impossible. 

In this paper we show that in fact QND measurements 
of the spin projection on the direction of the applied 
field B ]] z are possible with the use of spin dependent 
Josephson tunneling at low temperatures, when the ef- 
fect of quasiparticles is negligible. In this case the ampli- 
tude of Josephson oscillations depends on the state of the 
spin just before the measurement. The amplitude can be 
measured, e.g., in a circuit containing Josephson junc- 
tion with the spin, a dc-voltage source V and a resistor 
R, see Fig. 1. The dissipative spin-dependent current in 
the circuit appears when the voltage exceeds a threshold 
and, effectively, it measures the squared amplitude of the 
Josephson-current oscillations. 

The measurement can begin either when the voltage 
exceeding a threshold [12] or the tunneling in the pres- 
ence of such a voltage is switched on. The first possibility 
appears more natural for mesoscopic circuits, while the 
second could be realized with the use of the STM with 
a superconducting tip and a molecule with spin on the 
superconducting substrate. As the tip approaches the 
molecule, tunneling is turned on. 



FIG. 1. The circuit containing Josephson junction with a 
spin. 

The system we consider, i.e. a spin in a Josephson 
junction with applied voltage, is interesting also from the 
general viewpoint of quantum measurements. The entan- 
glement of the spin with the measurement apparatus is 
realized here by use of the non-dissipative macroscopic 
quantum system as an intermediate step, namely with 
the superconducting phase degree of freedom. Remark- 
ably, in such a system the amplitude of Josephson oscil- 
lations at frequency ujj — 2eV/h at finite bias V carries 
information on the initial state of the spin even for junc- 
tions made of singlet superconductors. Due to the spin 
conservation in a singlet-Cooper-pair tunneling in the ab- 
sence of quasiparticles (T 0) the average spin compo- 
nent, (Sz), is preserved to second order in the tunneling 
amplitude at y < 2Ao (Aq being the superconducting 
gap), though the spin operator Sz does not commute with 
the Hamiltonian 7i of the system. The noncommutativ- 
ity of Ti. and S enables the tunneling measurement of the 
z-component of the spin, while due to the preservation 
of (Sz) it is a QND measurement. 

We consider a Josephson junction where the spin- 
independent tunneling is described by amplitude To and 



the spin-dependent tunneling by Tg 
nian of the system is 



To. The Hamilto- 
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'H = Ho + 'Ht, Ho = 'Ha + 'Hb - jJ^BzSz, (1) 

a, a' 

+ Ts{cr) 

a, a' 



n,m,a,Q' 



where T^o and are the Hamiltonians of the supercon- 
ducting leads a and b, and aj^^ {ana) creates (annihilates) 
au electron in the lead a in the state n with the spin a. 
Further, a represents the three Pauli matrices, the local- 
ized spin 1/2 is described by the operator S, and /z is the 
magnetic moment of the spin. In the following we use 
notations B for iiB^ and lol = B/h. 

As we will show in the following, the Josephson cur- 
rent I{t) oscillates at the frequency ojj = 2eV/H, and 
the average oscillation amplitude depends on the state of 
the spin at t = 0, i.e. on |a|^ for the spin wave func- 
tion *s(0) = a| t) + /3| i)- We show also that the de- 
component of the current in the circuit with a resistor 
depends quadratically on the oscillation amplitude, and 
thus is sensitive to the value Sz at time t = 0. 

First, we derive the expression for the Josephson cur- 
rent across the junction as a functional of the phase dif- 
ference (f{t) across the junction and the spin wave func- 
tion 5's(0) at low temperatures, neglecting the effect of 
quasiparticles, which will be estimated later. We neglect 
quantum fluctuations of the phase difference across the 
junction, assuming that the charging energy (2e)^/C is 
lower than the Josephson energy ~ TqP^Aq, where C is 
the junction capacitance and p is the density of states 
per spin in the leads. We assume that T^p^.T'^p^ <^ 1 
and 2eV,B <^ Aq, and that all relevant frequencies are 
much lower than Aq / h. We find the Heisenberg operator 
for the current using the operator perturbation theory 
[13] with respect to the tunneling Hamiltonian Ht, and 
then we average this operator over the electron ground 
state of the Hamiltonian Ha + Hb- The current oper- 
ator {I{t))e obtained in this way is an operator in the 
spin space and ('i'*(0)|(/(t))e|*I's(0)) defines the average 
current for any initial state (0). Yet the measurement 
produces one of the two values, corresponding to the spin 
up/down states, Sz = ±1/2, rather than the average. 

At time t > the current is given by 



I{t) = exp(| m)7exp(-^ Ht) , 
exp(— — Ht) = exp(— — 7^oi)Texp 



nJo 



HT{T)dT 



(3) 
(4) 



1= -f-aL[ToSc,a'+T,{&)„,c.'-S]bmc,>+h.c. , (5) 



where we accounted for isotropic exchange interaction of 
tunneling electrons with localized spin, T is the time or- 
dering operator and A{t) — exp{iHot/h)Acxp{—iHot/h) 
is the operator in the interaction representation, A is the 
Schroedinger operator and A{t) is the Heisenberg oper- 
ator. After averaging over electronic degrees of freedom 



in the absence of quasiparticles and assuming a classical 
phase difference on the junction ip{t) wc find 



{m)e = ^ 



dt' sin 



fit) + fit') 



{F+{t - t')F{t - t')[n - TfSit) ■ S{t')] - h.c] , (6) 

where S(f) = exp(— 2CJiS'^t)S cxp(icij/^52f), and the 
Gor'kov Green functions are related to the Bessel func- 
tions as 

F+{t) = -F{t) = (7rAop/2)[Jo(Aot/^) - iNo{Aot/h)] . 

For singlet Cooper pairing the dependence of the current 
on the spin is isotropic due to the isotropic exchange cou- 
pling of the spin and the tunneling electrons. We see that 
the average current is expressed via the spin correlation 
function {**(0)|S(t)-S(t')|*s(0)). This function depends 
on the initial spin state, thus leading to the correspond- 
ing dependence for the supereurrent, because 

S{t) ■ S(0) = {2cosLVLt + 1)1/4: - iSzSmuLt, (7) 

where 1 is the unit matrix in the spin space. 
The functions F+{t)F{t), F+ {t)F{t) cos lolI and 
F~^ {t)F{t) sin u)Lt oscillate and drop on the scale Aq^ 
at _B ^ Aq. For the function <f{t) with characteristic 
frequencies well below Aq, we take (pit) « fit') and for 
t ^ Aq ^ we obtain 



(/(t)}e = /cSin^(t), ic = Ioi+IsSz. 



(8) 



-^(t^--t^]p^Ao, h 



4e, 

'h 



T^p^B. 



(j(t))e is the renormalized current via the Josephson 
junction obtained after integration over all high frequen- 
cies > Aq. It corresponds to the Josephson energy 

Ej{^{t),Sz}^ih/2e)(hi + IsSz)[l-cos^it)l (9) 

which can be considered as an effective Hamiltonian for 
the Josephson junction with the spin. The spin sup- 
presses the Josephson energy due to the negative con- 
tribution (3/4)Tg^ to /q which can lead to the formation 
of a TT-j unction when /q changes sign [7]. 

The current amplitude depends on the initial state of 
the spin with different values iov Sz = 1/2 and Sz = —1/2 
until the spin flips due to quasiparticles or due to re- 
laxation. Otherwise, on average, the spin does not flip 
in Josephson tunneling because the flrst electron of the 
Cooper pair may flip it, but then the second one restores 
the spin since the Cooper pair cannot carry spin; thus Sz 
on average is an integral of motion. 

The effect of the Josephson current on the spin can 
also be understood from the following argument. One 
sees from Eq. (2) that the spin feels the fluctuating ex- 
change field induced by tunneling electrons. 
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E 



aia{t){<^)a,a'bm.a'{t) + h.C. . (10) 



This field causes deviations of spin from the 2;-axis orien- 
tation due to presence of transverse components ax and 
ay when the first electron of a Cooper pair tunnells via 
spin. This spin response interferes with the spin depen- 
dent part of the current operator (5) giving rise to the 
spin dependent contribution in the current. 

Having obtained linear dependence of the Josephson 
critical current on the spin ^-component we have reduced 
the problem to the well studied case of a dc-SQUID mea- 
suring a flux qubit [4]. Here, for completeness, we provide 
a simplified version of the derivation of the dephasing and 
measurement times, valid for i? <C /i/e^ and T <C V. 

Without the resistor {R = 0) the system is dissipa- 
tionless, (fi{t) = ujjt and no measurement is performed. 
Introducing dissipation allows us to measure the ampli- 
tude of the Josephson-current oscillations, and, thus, the 
initial state of the spin. We consider the circuit in Fig. 1 
with R^O. The circuit dynamics is described by 



2e ^ 2eR 



Ic sin if 



(11) 



where the Langevin force ^(t) is the random current 
(Nyquist noise) generated by the resistor. Let us con- 
sider the simplest case when the junction capacitance 
can be neglected. This is possible if RCluj <g; 1. When 

V <C IcR no voltage drops on the junction and the dc- 
current is not sensitive to the value of the Josephson 
critical current and, thus, to the state of the spin. For 

V ^ IcR, solving Eq. (11) and averaging one obtains 
(fi{t) = u)jt + (IcR/V) cos{u)jt) (see e.g., [8]) and 

I = lcsm(p^ Ic sm{u! jt) + {llR/V) cos^ {ujjt) . (12) 

The dc-component of the current is given by 



f llR _ {Il + PJ4.)R . , JqV.^ 



nR . 



2V 



2V 



V 



(13) 



Now the dc-current is sensitive to the state of the spin. 
For a single measurement one of the two possible values of 
the current, (7o ± -fs/2)^-R/2y, is realized. In a multiple 
set of measurements these values occur with probabilities 
P+ = |a|2 and P_ = 1 - |a|2. 

Now we are in a position to derive the measurement 
rate Tm needed to resolve two values of the current corre- 
sponding to 5z = ±1/2 on the background of shot noise 
of Cooper pairs tunneling incoherently due to the resis- 
tor. The signal is 51 = IqIsR/V. The noise power of the 
background current Si{u)) is defined as 



1 r 

2 7-< 



dtexp{iwt){[I{t)I{0) + I{0)I{t)]) . (14) 



Si{ij) « 2eldc « el^R/V, a; < F < Aq, 



(15) 



assuming /q » /«. This expression can be obtained [8] 

from the Langevin equation (11) with the appropriate 
high frequency spectrum of £,{t) [9]. It describes the shot 
noise of individual Cooper pairs tunneling incoherently 
through the junction and dissipating each energy 2eV 
into the microscopic modes of the resistor [10]. We de- 
fine the measurement rate as: 



i6ir/{8Sj) « {I^R)/{8eV) 



(16) 



Note that the current correlation function for the shot 
noise is expressed via the phase correlation function as 



(mm) ^I^{sm^{t) sin cp{0)) . 



(17) 



Next we discuss the back-action effect of the measure- 
ment, i.e. the effect of incoherent Cooper pair tunnel- 
ing, on the spin dynamics. It leads to the dephasing of 
transverse spin components, i.e. to the destruction of 
coherent superposition of states Sz ~ ±1/2. To derive 
the dephasing rate, F^, we consider the x-component of 
spin averaged over the electronic degrees of freedom. We 
obtain 

h ~ /■* 

{Sx{t))e = rjit)Sx{t) + Sy(t) J COSip{T)dT , (18) 

where ri{t) < 1 is an oscillating function close to unity 
[see Eq. (20)], describing the reduction of the spin 
amplitude. Its origin will be discussed later. Here 
we focus on the second term in the right hand side. 
This term describes the phase accumulation due to 
the effect of a random magnetic field acting on spin, 
hzr{t) = [hls/{2e)] cos(p{t), in agreement with the ef- 
fective Hamiltonian (9). The dephasing rate is given by 



^ _J_ fhis 

2^12 \ 2e 



dr {cos (p{t) cos (f {0)) . (19) 



At y > IcR it is given by 



By use of Eqs. (15), (17) we get F^ « P^R/{8eV) « F„. 

The inequality F^ > F^ is the fundamental property 
of quantum mechanics as it follows from the uncertainty 
relations (see, e.g., review [11] and references therein). It 
was proven for the case of a linear detector with the di- 
rect coupling of the detector to the measured observable, 
see e.g., Ref. [4,14]. In our setup the interaction of spin 
with tunneling electrons is nonzero only in the second 
order of perturbation theory. However, effectively we 
obtain a linear relation between the spin ^-component 
and the superconducting current accessible experimen- 
tally, see Eqs. (8), (9). Consequently, the optimal rela- 
tion Frf w F„ is satisfied for our measurements. 

Now we discuss the effects which are beyond the accu- 
racy of the effective Hamiltonian (9). We calculate the 
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z-component of the spin, {Sz{t))e averaged over the elec- 
tronic degrees of freedom, as was done in Eqs. (3), (4). 
We derive for t^h/Ao: 

(4(i))e = (exp(im)5,exp(-zHf))e = vit)Sz, 
r]{t) = l- 2T>2 [1 + cos ^{t)] . (20) 

The maximum average value of the spin 2;-component 
is slightly smaller than 1/2. This is because spin 
flips can happen after one electron of a Cooper pair 
passes the spin. The time of passing is ~ h/Ao and 
probability to scatter on the spin per unit time is ^ 

p'^ {Ao /h) cos ip{t). Also the second term in Eq. (18), 
which corresponds to an extra phase in the spin rotation, 
can be understood in terms of the spin flips. Each pas- 
sage of a Cooper pair changes this phase by Aa ~ B/Aq 
and in the case of coherent tunneling during the time 
t this phase changes by a{t) ~ T"^ {Aq / h){B / Ao)t ~ 
Tgp^LOLt- Hence, in addition to B there is the efl'ec- 
tive magnetic field acting on spin hz ~ T^p^B as follows 
also from the effective Hamiltonian (9). For uncorrelated 
tunneling of Cooper pairs this field is random and in ad- 
dition to regular change it leads to phase diffusion as 
{[a{t) - (a(t))]a(0)) oc {Aa)H oc BH and oc B^. 

For comparison, we consider a normal tunnel junction 
with unpolarized electrons in the leads. We find the cur- 
rent sensitivity to the spin z component also, 

Ottp 

{Iit)),^—p'[{2Ti + -iT])Vl- (21) 

T^{\V + B\-\V -B\)Sz]. 

Now, however, the spin flips are allowed because tun- 
neling electrons are not paired. For short times 
h[T'^plV]-'^ > t > nV-'^ we obtain 

{Sz{t))e = Sz- 27rT2p'[24max(|y|, \B\) - B]t, (22) 

which gives at eV ^ B the spin-flip (demolition) rate 
r|j ~ TgP^eV/h. The dcphasing rate is about the same. 
The measurement rate given by HTm ~ p'^ B'^ / {Tq eV) 
is smaller than the spin-flip rate. Hence, tunneling mea- 
surement of spin is impossible in the normal state. 

Thermal quasiparticles cause spin flips as the electrons 
do in the normal state. The demolition rate is propor- 
tional to the fraction of quasiparticles. At T ^ Aq it is 
exponentially smaU, hVl < T^p'^eV e^p{-Ao/T). 

In this paper we have considered spin 5* = 1/2. We 
note that also for arbitrary S the current operator {I{t))e 
is a linear combination of operators S'", < n < 2S. 
It depends on the spin initial state and differs for the 
states (Sz) — ±5" in contrast to the classical spin (dis- 
cussed in Rcf. [3]) for which it depends only on S"^ be- 
cause S{t) ■ S(0) = (5*2 - S^z) cosLULt + si 

We note, that other strategies of spin measurement 
could also be considered. E.g., a Josephson junction with 
a spin could be inserted into a superconducting ring, a 



voltage across the junction could be induced by chang- 
ing magnetic flux across the ring and the alternating flux 
induced by Josephson oscillations could be measured by 
a secondary circuit coupled inductively to the supercon- 
ducting ring. This will be studied elsewhere. 

In conclusion, we have shown that the nondemolition 
measurements of a quantum spin are possible using the 
Josephson oscillations. The average amplitude of Joseph- 
son oscillations depends on the initial state of the spin af- 
ter the tunneling or voltage is switched on and the voltage 
exceeds the dissipative threshold. This amplitude can be 
found measuring the dc-current in the circuit with a resis- 
tor. Use of the singlet-pair tunneling as an intermediate 
nondissipative system allows us to obtain information on 
the initial value of the z component of the spin without 
flipping the spin. Thus a QND tunneling measurement 
of spin is performed. 
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